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Abstract 

Mathematical  development  of  preliminary  performance/design  equations  for  a  hydrogen/air,  solid  polymer  acid  electrolyte  based  fuel  cell  is 
presented.  The  development  is  based  on  the  principles  of  transport  phenomena,  intrinsic  electrochemical  kinetics,  and  classical  thermodynamics. 
The  developed  formulation  is  intended  to  quantitatively  describe  the  mass  fraction  profiles  of  the  chemical  species,  hydrogen  and  oxygen,  in  the  cell 
anode  and  cathode  diffusion  and  electrocatalytic  reaction  layers  as  a  function  of  the  distance  in  the  proton  transport  direction  at  an  axial  distance 
parallel  to  the  cell  anode  or  cathode  channel  flow.  Given  the  cell  geometry,  chemical  species  and  charge  transport,  and  intrinsic  electrochemical 
kinetic  parameters,  the  developed  formulation  can  be  employed  to  compute  the  species  local  mass  fluxes  and  predict  the  cell  anode  and  cathode 
cell  overvoltages  for  a  desired  geometric  current  density.  The  presented  single  cell  performance  predictive  formulation  has  also  been  linked  to  the 
formulation  needed  to  predict  the  performance  of  a  stack  of  a  number  of  identical  PEMFCs  connected  in  series. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  current  need,  throughout  our  world,  is  power  generation  with  greater  energy  efficiency  and  environmental  protection.  In 
general,  fuel  cells  are  a  clean,  efficient  and  quiet  source  of  power  generation.  The  polymer  electrolyte  membrane  fuel  cell,  with 
hydrogen  as  the  primary  fuel,  continues  to  be  developed  due  to  its  low-temperature  operation,  rapid  start-up  and  high  power 
density.  The  polymer  electrolyte  membrane  fuel  cell  (PEMFC;  also  called  the  proton  exchange  membrane  fuel  cell)  converts 
the  chemical  energy  of  hydrogen  and  oxygen  from  air  directly  into  electrical  energy,  with  heat  and  water  as  by-products.  In  the 
near  future,  the  PEMFCs,  using  hydrogen  and  air  as  reactants,  are  expected  to  replace  the  less  efficient,  pollutant  (e.g.,  CO,  NOv, 
unburned  hydrocarbons,  soot)  generating  internal  combustion  engines  burning  fossil  fuels  acquired  from  dwindling  terrestrial  energy 
sources.  Currently,  the  PEMFCs  are  not  a  commercially  viable  alternative  to  the  internal  combustion  engines.  However,  an  enhanced 
understanding  of  the  various  transport  and  electrochemical  processes  that  occur  within  a  PEM  fuel  cell,  through  mathematical 
modeling  and  experimental  work,  can  help  in  the  acceleration  of  commercialization  of  the  PEMFC  technology  by  providing  gains 
in  cost  reduction,  energy  conversion  efficiency,  and  power  output. 

Recent  one-dimensional  models  of  the  PEMFC  are  given  in  Refs.  [1-8].  These  models  provide  useful  insight  and  reasonable 
predictions  of  the  cell  performance  in  the  low  and  intermediate  current  density  ranges,  but  are  not  able  to  simulate  the  rapid  electric 
potential  drop  observed  experimentally  at  higher  current  densities.  Two-dimensional  models  are  given  in  Refs.  [9-21].  Recent  three- 
dimensional  works  are  found  in  Refs.  [22-33].  Most  of  these  models  found  in  the  literature  are  shown  to  simulate  the  overall  fuel 
cell  performance  (i.e.,  cell  voltage  versus  current  density)  reasonably  well.  This  agreement,  irrespective  of  the  dimensionality  of  a 
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Nomenclature 

aa  activity  of  species  a 

ae ff  platinum  surface  area  fraction  effective  for  electrochemical  reaction 
aA  active  electrocatalyst  surface  area  per  unit  ARL  volume  (m-1) 

ac  active  electrocatalyst  surface  area  per  unit  CRL  volume  (m-1) 

Aj,  Bj,  Cj,  Dj  heat  capacity  coefficients  for  an  ideal  gas  i 
c  total  gas  phase  molar  concentration  (mol  m-3) 

ca  gas  phase  molar  concentration  of  species  a  (mol  m-3) 

Da  gas  phase  mass  diffusivity  of  species  a  (m2  s-1) 

E  cell  voltage  (V) 

Eq  reaction  activation  energy  (J  mol-1) 

/i  defined  by  Eq.  (39) 

F  Faraday’s  constant  (96,487  C  g  equiv.-1) 

FA  defined  by  Eq.  (31b) 

Ff  defined  by  Eq.  (82b) 

AG  Gibbs  free  energy  change  (J  mol- 1 ) 

AH  enthalpy  of  reaction  (J  mol- 1 ) 

AHj  enthalpy  of  reaction  given  by  Eq.  (3)  at  temperature  T  (J  mol-1) 

A  Ht  thermal  energy  production  by  Eq.  (3)  for  combustion  (W) 

;geom  geometric  current  density  (A  m-2) 

z'o.s  intrinsic  exchange  current  density  (Am-2) 

It  total  current  (A) 

kc  mass  transfer  coefficient  (m  s-1) 

kc  electronic  conductivity  (G- 1  m-1) 

kb  frequency  factor  for  hydrogen  electrochemical  reaction  at  the  electrocatalyst  surface  (mol  s-1  m-2  N  m-2) 

L  length  (m) 

mpt  geom  platinum  loading  (kg  m-2) 

Ma  molecular  weight  of  species  a  (kg  mol- 1 ) 
n  number  of  cells  in  a  series  connected  stack 

ha  mass  flux  of  species  a  (kg  m-2  s- 1 ) 

Na  molar  flux  of  species  a  (mol  m-2  s-1) 

Ain, «  total  molar  flow  rate  in  of  species  a  (mol  s-1) 

Aexit.a  total  molar  flow  rate  out  of  species  a  (mol  s-1) 

Aexit.t  total  molar  flow  rate  out  (mol  s-1 ) 

—  A  ;V||2  hydrogen  consumption  rate  ( mol  s- 1 ) 

Q remove  thermal  energy  needed  to  be  removed  to  remain  isothermal  (W) 
r pt  platinum  average  spherical  particle  radius  (m) 

ra  electrochemical  molar  reaction  rate  of  species  a  per  unit  of  catalyst  surface  (mol  s- 1  m-2) 

Ra  electrochemical  molar  reaction  rate  of  species  a  per  unit  geometric  volume  (mol  s-1  m-3) 

Ra, mass  electrochemical  mass  reaction  rate  of  species  a  per  unit  geometric  volume  (kg  s-1  m-3) 
s  platinum  specific  surface  area  (m2  kg-1) 

pa  partial  pressure  of  species  a  (N  m-2) 

P  gas  mixture  total  pressure  (N  m-2) 

PEMFC  polymer  electrolyte  membrane  fuel  cell  or  proton  exchange  membrane  fuel  cell 
R  universal  gas  constant  (8.314  J  mol-1  K-1  or  8.314Nm-2  m3  mol-1  K-1) 

T  Temperature  (K) 

W  width  (m) 

^electric  electric  power  (W) 

x  coordinate  axis  for  cell  length  or  distance  along  x  axis  (m) 

y  coordinate  axis  for  cell  width  or  distance  along  y  axis  (m) 

ya  mole  fraction  of  species  a 

Z  coordinate  axis  for  cell  thickness  or  distance  along  z  axis  (m) 

^near-side  5ACH  +  SADL  (m) 

4r-ride  5ACH  +  SADL  +  SA  (m) 
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,CRL 

^-ns 

,CRL 

>L 

3ACH  +  ^ADL  +  ^A  +  5SEP(m) 

3ACH  +  ^ADL  _j_  ^A  _|_  ^SEP  +  (m) 

3ACH  +  ^ADL  +  SA  +  ^SEP  +  (m) 

VCDL 

zfs 

3ACH  +  ^ADL  +  SA  +  ^SEP  +  &C  +  5CDL(m) 

Greek  letters 

a 

defined  by  Eq.  (35c) 

ai 

defined  by  Eq.  (84b) 

aa 

anodic  charge  transfer  coefficient 

<*c 

cathodic  charge  transfer  coefficient 

8 

thickness  (m) 

£ 

volume  fraction  or  void  fraction 

K 

defined  by  Eq.  (35b) 

£h2o 

moles  of  water  transported  to  CRL  per  mole  of  H+  migrating  from  ARL  to  CRL 

t] 

voltage  loss  (V) 

^thermal 

thermal  efficiency 

Vi 

stoichiometric  coefficient 

p 

gas  phase  density  (kg  m-3) 

O 

ionic  conductivity  (ohm-1  m-1) 

rp(g) 

tortuosity  factor  of  gas-filled  pores 

<p 

electric  potential  in  ARL  (V) 

A  tp 

Galvanic  electric  potential  difference  in  ARL  (V) 

(P 

electric  potential  in  CRL  (V) 

A  (p 

Galvanic  electric  potential  difference  in  CRL  (V) 

COot 

mass  fraction  of  species  a 

Superscripts 

A 

anode  reaction  layer 

ACH 

anode  channel 

ACH-DL  anode  channel  +  anode  diffusion  layer 

ACH-DL-RL  anode  channel  +  anode  diffusion  layer  +  anode  reaction  layer 

ADL 

anode  diffusion  layer 

ARL 

anode  reaction  layer 

cell 

cell 

contact 

contact  resistance 

C 

cathode  reaction  layer 

CCH 

cathode  channel 

CDL 

cathode  diffusion  layer 

CRL 

cathode  reaction  layer 

m 

reaction  order 

sat’d  vapor  saturated  vapor 

stack 

number  of  n  cells  in  series 

SEP 

solid  polymer  electrolyte  membrane 

t 

total 

O 

standard  state 

Subscripts 

ave 

average 

carbon 

carbon  material  used  in  electrodes 

e“ 

electronic 

eff 

effective 

electrolyte  solid  polymer  electrolyte 

fs 

far  side 

g 

gas 

geom 

geometric 
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H+ 

protonic  or  ionic 

ns 

near  side 

rev 

reversible 

PU-,,0 

reference  partial  pressure  of  hydrogen 

Pt 

bulk  platinum  catalyst 

t 

total 

z 

coordinate  axis  for  cell  thickness 

T 

cell  temperature 

To 

reference  temperature  (298.15  K) 

0 

reference  condition 

a 

H2,  N2,  CO2,  H2O,  etc. 

a-SACU 

chemical  species  a  at  the  ACH/ADL  interface 

P 

H2,  N2,  CO2,  H2O,  etc. 

Q 

ohmic 

model,  is  quite  likely  due  to  the  adjustment  of  one  or  more  parameter  values.  A  three-dimensional  modeling  work  [34]  has  been 
validated  against  local  current  density  data  obtained  experimentally  with  a  segmented  cell.  All  these  models  attempt  to  model  a 
PEMFC  in  a  realistic  fashion  without  resulting  in  too  much  of  computational  expense. 

A  mathematical  formulation  is  developed  below  for  a  hydrogen/air  PEMFC  utilizing,  for  example,  Nafion,  PBO(poly(p- 
phenylene-2,6-benzobisoxazole))/PBI  (poly (benzimidazole))  based  phosphoric  acid,  or  sulfonated  polyarylene  ether  or 
Nafion-Teflon-Zr(HP04)2  as  an  electrolyte.  This  two-dimensional  model,  originally  developed  for  a  PEMFC  fueled  by  hydro¬ 
gen  gas,  can  be  extended  to  predict  the  performance  of  PEMFC  fueled  by  a  gas  mixture  containing  hydrogen,  carbon  monoxide, 
small  amount  of  oxygen,  carbon  dioxide  and  nitrogen  as  found  in  the  gas  mixture  issuing  from  a  reformer  dedicated  for  the  production 
of  fuels  for  fuel  cell  applications  [35].  Isothermal,  steady-state  conditions  have  been  assumed. 

2.  Mathematical  formulation 

Fig.  1  shows  the  schematic  of  a  single  proton  exchange  membrane  fuel  cell  (PEMFC).  Operation  of  such  a  PEMFC  is  assumed  to 
be  under  the  isothermal  and  steady-state  conditions.  Hydrogen  gas  or  gas  mixture,  containing  hydrogen,  nitrogen,  carbon  dioxide, 
etc.,  from  a  reformer  enters  the  anode  channel  and  flows  in  the  x-direction.  At  any  x-location,  hydrogen  diffuses  from  the  anode 
channel  bulk  flow  to  the  porous  carbon  anode  diffusion  layer  (ADL).  Then,  it  diffuses  through  the  ADL  to  the  interface  between  the 


Fig.  1.  Schematic  of  a  single  fuel  cell.  Note:  The  x-direction  is  J_  the  y-z  plane.  Fuel  and  oxidant  gas  streams  in  the  anode  and  cathode  channel  are  parallel  to  the 
x-direction.  LACH,  LCCH:  anode  and  cathode  channel  length,  respectively,  parallel  to  the  x-axis.  <5ACH,  5ADL,  <$A,  <$SEP,  8C,  <5CDL,  and  <$CCH  are  the  dimensions  of 
anode  channel,  anode  diffusion  layer,  anode  reaction  layer,  solid  polymer  electrolyte  membrane,  cathode  reaction  layer,  cathode  diffusion  layer,  and  cathode  channel, 
respectively,  parallel  to  the  z-axis.  WACH  and  WCCH  are  the  anode  and  cathode  channel  dimensions,  respectively,  parallel  to  the  y-axis. 
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ADL  and  the  anode  reaction  layer  (ARL).  The  porous  reaction  layer  is  assumed  to  be  composed  of  the  catalytic  platinum  particles, 
in  contact  with  an  acid  solid  polymer  electrolyte  such  as  Nafion  or  PBO/PBI  electrolyte,  dispersed  on  the  porous  carbon  granules. 
Hydrogen  diffuses  toward  the  interface  between  the  ARL  and  solid  polymer  electrolyte  separator  layer  sandwiched  between  the  anode 
and  cathode  reaction  layers.  In  the  ARL,  hydrogen,  while  it  diffuses,  is  consumed  via  the  overall  electrochemical  reaction  given  by 
Eq.  (1).  The  generated  electrons  enter  an  external  electric  circuit  to  pass  through  an  external  electrical  load.  The  generated  protons, 
H+,  migrate  internally  through  the  electrolyte  separator  to  the  porous  cathode  reaction  layer  (CRL)  where  they  cause  reduction  of 
oxygen  in  the  presence  of  electrons  that  return  from  the  external  electric  load  to  the  cathode  reaction  layer. 

Atmospheric  air,  containing  varying  degrees  of  water  vapor,  enters  the  cathode-side  flow  channel.  At  any  x-location  in  the  flow 
channel,  oxygen  diffuses  toward  the  cathode  reaction  layer  in  a  negative  ^-direction.  Oxygen  first  diffuses  through  the  chemical 
species  concentration  boundary  layer  to  the  porous  carbon  cathode  diffusion  layer  (CDL).  The  chemical  species  concentration 
boundary  layer  may  or  may  not  prevail  at  the  porous  surface  of  the  CDL  depending  on  the  oxygen  reduction  rate  in  the  CRL 
and  the  cathode  channel  flow  Reynolds  number  value.  It  then  diffuses  through  the  CDL  to  the  interface  between  the  CDL  and 
relatively  thinner  CRL.  The  porous  CRL  is  assumed  to  be  composed  of  the  catalytic  platinum  particles,  in  contact  with  an  acid 
solid  polymer  electrolyte,  for  example,  Nafion  or  PBO/PBI,  dispersed  on  the  porous  carbon  granules.  In  the  CRL,  oxygen  diffuses 
toward  the  interface  between  the  CRL  and  solid  polymer  acid  electrolyte  separator  layer,  while  it  is  simultaneously  consumed  via  the 
following  overall  electrochemical  reduction  reaction  given  by  Eq.  (2).  The  electrochemically  produced  water  in  the  cathode  reaction 
layer  diffuses  through  the  cathode  reaction  and  diffusion  layers.  It,  then,  diffuses  through  the  species  concentration  boundary  layer 
prevailing  at  the  porous  CDL  surface  to  join  the  cathode  channel  bulk  flow.  It  is  here  noted  that  the  oxygen  mass  transport  is  in  the 
direction  of  decreasing  z  while  water  vapor  transport  is  in  the  opposite  direction.  Transport  of  the  gaseous  species  through  the  solid 
polymer  electrolyte  separator  layer  is  assumed  to  be  negligibly  small. 

Given  below  are  the  main  mathematical  model  equations.  One  may  find  the  detail  behind  the  formulation  development  elsewhere 
[36], 

The  half  cell  reactions  at  the  anode  and  cathode  electrocatalytic  reaction  layers,  respectively,  are: 

H2^  2H+  +  2e“  (1) 

2 02(gj  +  2H+  +  2e“  H20(g)  (2) 

The  overall  cell  reaction  is 


H2(g)  +  ^02(g)  H20 


The  reversible  cell  voltage  is  given  by 


c-cell 

rev 


=  E°  + 


RT 
2 F 


In 


(Pn2/p°)(po2/p°)l/2 

«h2o 


(3) 

(4) 


In  the  denominator  of  Eq.  (4),  the  constant  two  corresponds  to  2  mol  of  H+  or  e  per  1  mol  of  H2  consumed  electrochemically.  The 
standard-state  reversible  cell  voltage,  E°,  at  the  cell  temperature,  T,  is  given  by 


E°  = 


-  AGj 
2  F 


(5) 


For  convenience,  A Gj,  can  be  computed,  with  negligibly  small  error,  from  the  following  equation  since  it  can  be  assumed  that 
(3AG  j/dT)v  is  approximately  constant  over  the  range  of  interest  [37]: 


AG°r  =  (  ^  (A  G°Tn  -  A  HZ)  +  A  H\ 


To) 


Do 


(6) 


Total  ohmic  voltage  drop  in  the  electrolyte  separator  and  anode  and  cathode  diffusion  and  electrocatalytic  reaction  layers  is  given  by 

4  =  4EP  +  (4DL  +  4DL)  +  (4  +  4)  +  7contact  (7) 

The  ohmic  voltage  loss  due  to  the  resistance  to  H+  transport  in  the  electrolyte  separator  (SEP)  layer,  ;;^p  is  computed  from: 

<5sep 


„SEP  -  7 
—  1 


geom  " 


crH+  ,e£f, electrolyte 

The  ohmic  voltage  drop  in  the  anode  and  cathode  diffusion  layers  are  given  by 
5adl 


(8) 


ADL 


‘geom  ADL 

'V.eff 


(9a) 
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where 


lAUL  _  AUL  i.AUL 

Ke-,eff  fccarbon'''e_,eff,carbon 


(9b) 


and 


xCDL 

CDL  _  ■  _ 

rin  -  Igeom  CDL 

fj-.eff 


(10a) 


where 


.CDL  _  CDL  .CDL 
A'e_,eff  fccaibon'ce_,eff, carbon 

The  ohmic  drop  in  the  anode  reaction  layer  due  to  the  resistance  to  H+  and  electron  conduction  is  given  by 
„A  _  A  4-  nA 


(10b) 


(11a) 


igeom^  /  1 


9  \  nA  uA 

V  H+,e£f  'V.eff, 


(lib) 


where 


A  _  A  A 

H+  ,eff  electrolyte  aH+ ,eff, electrolyte 


and 


(11c) 


eff 


=  e; 


A 


/.A 


carbon^e-  ,eff, carbon 


(lid) 


The  factor  2  in  the  denominator  of  the  right-hand  side  expression  of  the  equal  sign  in  Eq.  (lib)  appears  to  account  for  the  fact  that 
the  H+  flux  is  zero  at  z  —  ^near-side  just  on  die  ‘minus  side’  of  the  interface  between  the  anode  diffusion  and  reaction  layers  and 
it  corresponds  to  /ge0m  at  z  =  z^r  Siic|e,  just  on  the  ‘plus  side’  of  the  interface  between  the  anode  reaction  and  polymer  electrolyte 
separator  layers.  The  electron  flux  is  zero  at  z  =  (zj'lr_SK|c)  and  it  corresponds  to  /gcom  at  z  =  (^near-side)  •  This  means  that  the 
arithmetic  average  current  density,  igeom/2,  has  been  used  in  Eq.  (lib)  to  estimate  r fy. 

The  ohmic  voltage  drop  in  the  cathode  reaction  layer  due  to  the  resistance  to  H+  and  electron  conduction  is  given  by 


c  c  ,  c 
In  —  ^h+,12 


(12a) 


4geom 


5C 


+ 


H+.eff 


kC 

e-.e£f/ 


(12b) 


where 


c  _  c  c 

^H+.eff  ^electrolyte^H+.eff,  electrolyte 


(12c) 


and 


e_,eff  ‘ 


c  carbon  e~  ,eff, carbon 


(12d) 


Voltage  losses  associated  with  the  contact  resistances  between  the  various  cell  layers  is  given  by  )/contact  and  is  due  to  the  lack  of 
‘perfect  contact’  and  is  a  function  of  stack  pressure.  Perfect  contact,  as  a  first  order,  is  assumed  in  this  analysis. 

2.1.  The  fuel  cell  anocle-side  formulation 

The  flux  of  a  chemical  species  a  in  the  z-direction  at  an  ^-location  through  the  fluid  phase  concentration  boundary  layer  prevailing 
at  the  porous  anode  carbon  diffusion  layer  (ADL),  in  general,  is  given  by 


iT  .ACH  ADL/  ACH  „  s 

Na,eff,z  ~  kc  a  £g  (Ca  -  Ca_s ACH  ) 


(13a) 


where  a  =  Hz,  N2,  CO2,  H2O,  etc.  The  mass  flux  of  a  chemical  species  a  in  the  z-direction  through  the  concentration  boundary  layer 
is  given  by 


ha.ciT.-  —  Na  e{{  z  Ma  —  k 


ACH  „ADL  j^ACH  „ACH,-,.,ACH 


c,a  ^ g 


K 


Ua-SAcil ) 


(13b) 
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(13c) 


.,ACH 


,  ,ACH  _ 
I = 


w  ACI I 


N: 


ACH 


v  Nach 

L^jol  ol 

ACH  */t 
ya  ma 

MACH 

pACH 


RT 


(13d) 

(13e) 

(13f) 


The  value  of  the  mass  transfer  coefficient,  kA™,  can  be  determined  using  standard  semiempirical  equations  [36,38]. 

The  gas  mixture  density  in  the  anode  channel  can  be  computed  from: 
pACH 

p  =  (14> 

OL 

The  mass  flux  equation  for  transport  of  a  chemical  species  through  the  porous  carbon  anode  diffusion  layer  is  given  as 
da>a 
dz 


na,zfi,z  —  P^a,eif  3" 

P 


(15a) 


where  a,  /?  =  H2,  N2,  CO2,  H2O,  etc.,  and 

nADL  .ADI. 


nADL  _ 
Uoifi{f  ~ 


'a- gas°g 
tADL 
TP(g) 


(15b) 


Correction  for  the  Knudsen  diffusion  effect  on  the  species  transport  in  the  pores  should  be  included  if  the  average  pore  radius  is  less 
than  lOnm  or  100  A. 

For  the  assumption  of  np.eff.z  —  0  (when  /3  =  N2,  CO2,  H2O),  with  a  =  FE,  Eq.  (15a)  reduces  to 


«H2,eff,z 


p/)ADL 

^H2,eff 

(1  —  &>H->) 


-dWa 

dz 


(16) 


Eq.  (16)  is  solved  for  t»\  1, ,  with  the  boundary  condition,  at  z  =  5ACH,  to h2  =  o)u  _t  sach  =  mass  fraction  of  hydrogen  in  the  gas  mixture 
at  the  porous  surface  of  the  carbon  ADL  facing  the  anode  channel  flow.  The  result  is 


o>h2  =  1  —  (1  —  <yH,il5ACH)exp 


»H2,eff,z 


(z-5ach) 


valid  for  5ACH  <  z  <  (5ACH  +  <5ADL).  Eq.  (13b)  is  used  to  eliminate  &>Hn  3ach  from  Eq.  (17)  to  obtain: 


ft>H2  =  1  —  <  1  —  4™ 


,ach£adlmachcach  I 


«H2,eff,i 


exp 


»H2,eff,; 

,P^H2Deff, 


(z  —  <5ach) 


(17) 


(18a) 


valid  for  <5ACH  <  z  <  (<5ACH  +  <5 ADL).  Eq.  (18a)  can  be  employed  to  calculate  the  mass  fraction  of  hydrogen,  to h9  ,  as  a  function  of 
distance  z  in  the  porous  carbon  anode  diffusion  layer  provided  «H2,eff,z,  corresponding  to  a  current  density,  ige om,  and  the  anode 
channel  bulk  gas  concentration  of  FE,  at  an  v-location  are  known. 

For  the  situation  of  non- zero  mass  flux  of  water  vapor  in  the  anode  diffusion  layer,  the  following  derived  expression  should  be 
used  to  compute  the  hydrogen  mass  fraction  profile  in  the  anode  diffusion  layer  [40]: 


&>h2  =  &  ■ 


fi)  _  .  ,ACH  , _ nH2,eff,z _ 1 

H2  kASHeADLMACHcACH  J 


exp 


ftH2,eff ,z 

KPD 


<9(z  -  «5ACH) 


(18b) 


where  ©  =  1/(1  +  18£h2o)- 

In  the  anode  reaction  layer  (ARL),  hydrogen  is  transported  as  well  as  consumed  via  the  electrochemical  reaction  at  the  active 
catalyst  platinum-electrolyte  interface  to  produce  H+  ions  and  electrons.  The  hydrogen  consumption  rate: 


(-4)  = 


‘O.s 

2  F 


exp 


<444 

RT 


(19) 
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where  the  anodic  intrinsic  exchange  current  density  is  given  by 


,-A 


•A 


O.s  <0.s,T0,Ph2,o  eXP 


zSL 
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1 
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f  P»2 


\Ph2,o 


(20) 


and 


1*7  I  —  T  Ercy  —  A(pi£Y\  —  |(<PPt  ^electrolyte)  (^Pt  ^electrolyte  )rcv 


(21) 


Note  that  when  an  impurity,  for  example  CO,  is  present  in  the  hydrogen  fuel  feed  mixture  above  10 ppm  [39],  it  acts  as  the  poison 
for  the  hydrogen  oxidation  or  de-electronation  reaction  on  the  cell  anode  surface.  Its  effect  to  decrease  the  rate  of  hydrogen  de¬ 
electronation  on  the  anode  surface  should  be  accounted  for  using  an  experimentally  determined  kinetic  deactivation  function  in  Eq. 
(19).  For  the  consumption  of  hydrogen  by  the  first-order  reaction,  m  =  1  in  Eq.  (20). 

If  it  is  assumed  that  the  catalyst  crystallites  are  spherical  and  their  average  size  is  i  pt  (radius),  the  platinum  surface  area  per  unit 
mass  of  platinum  crystallites  is  given  by 


4?r  (r£)“ 


(4/3)7r(r&)  pP,  'ptPPt 


(22) 


If  some  fraction  of  the  catalyst  surface  is  not  in  contact  with  the  electroactive  regions  of  the  polymer  electrolyte;  then,  the  reaction 
effective  electrocatalyst  interfacial  area  is  given  by 


4t : 


(23) 


For  a  platinum  loading  in  the  anode  reaction  layer  of  ;«pt,  the  electrochemical  reaction  active  electrocatalyst  surface  area  per  unit 
anode  reaction  layer  volume  is  given  by 


aA  —  4f»pt 


(24) 


The  relation  of  mg.  to  the  catalyst  loading  per  unit  geometric  area,  m 


Pt,geom 


,  of  the  anode  reaction  layer  of  thickness,  <5  ,  is  given  by 


,nPt,geom  ^  mPt 


(25) 


The  molar  hydrogen  consumption  rate  per  unit  anode  reaction  layer  volume,  — /?jj  ,  's  given  by 


-<  = 


*0,s 

2  F 


A 

a  exp 


ro^|?7A|F 


RT 


(26a) 


The  hydrogen  consumption  rate  in  mass  units  is  given  by 


-R 


0,s 


H2,mass 


ss  =  (-Rh2)Mh2  =  I  —  ]  «  Mu2  exp 


2  F 


RT 


(26b) 


For  the  assumption  of  hydrogen  transport  in  the  relatively  thin  anode  reaction  layer  via  molecular  mass  diffusion,  the  hydrogen  mass 
transport  flux  in  the  "-direction  in  the  anode  reaction  layer  is  given  by 


•  AR.L  _  rvARL 

nH2,eff,z  —  PuH2,eff 


— dtt>H2 

ds 


(27) 


A  mass  balance  for  hydrogen  over  a  thin  spatial  element  in  the  anode  reaction  layer  was  applied  to  obtain  the  following  differential 
equation  accounting  for  the  simultaneous  transport  and  consumption  of  hydrogen  via  electrochemical  reaction: 


fl.'.ARL 

Q”H2,eff.; 

dz 


‘0,s 


a  mh2  exp 


raA\rjA\F 


RT 


(28) 


On  substitution  for  n^fet r»  from  Eq-  (27)  into  Eq.  (28),  for  the  assumption  of  constant  (p^HTeff)  [40],  the  following  differential 
equation  was  obtained: 
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Using: 


_  ,,A  1-Eq 
2 F  0  £XP  (  RT 


Pn2 


=  I(q  exp 


0} h2 


Eq.  (29)  was  transformed  to  the  following  form: 
d2&>n2 


d  z2 


rA 

F\  ®h2 


where 


F?  = 


M 


ACHpACHflA 
n  7)ARL 


k[)  exp  - 


RT 


(30a) 

(30b) 


(31a) 


(31b) 


where  k ^  =  frequency  factor  for  the  hydrogen  de-electronation  or  oxidation  reaction  at  the  electrocatalyst  surface  for  the  assumed 
first-order  reaction.  Eq.  (31a)  was  solved  using  the  boundary  conditions:  at  z  =  <5ACH  +  <5ADL,  u> h,  =  3ach-adl>  the  hydrogen  mass 
fraction  on  positive  side  of  the  interface  between  the  anode  diffusion  and  reaction  layers;  atz  =  <5ACH  +  <5ADL  +  <5ARL,  (da)|i2/d")_  =  0. 
The  resulting  solution  to  compute  &>h2  as  a  function  of  the  distance  z  is  given  as 


(32) 


or,  further  simplification  yields: 


cou2  =  0) 


+ 

H2iaach-dl 


Eq.  (32)  is  valid  for  (5ACH_DL  <z  <  ^ach-dl-RL  [  so]ution  given  in  Eq.  (32)  was  obtained  using  the  assumption  that  variation  in 
|  r]A  |  with  respect  to  z  in  a  thin  reaction  layer  is  negligibly  small.  In  general,  activity  of  hydrogen  on  the  negative  and  positive  sides  of 
the  interface  between  the  anode  diffusion  and  reaction  layers  can  be  assumed  to  be  the  same  if  hydrogen  gas  is  present  in  different 
phases.  If  it  is  assumed  that  hydrogen  is  present  in  the  gas  phase  across  the  interface  between  the  anode  diffusion  and  reaction  layers; 
then,  it  is  reasonable  to  assume: 


ft,H2^ACH-DL  —  ®H2^ACH-DL 


(33) 


Using  Eq.  (18a),  the  following  equation  was  obtained: 


<wH2,,5ach-dl  —  ®H2,<5ach~dl 


_  ,  J  |  ,  ACH  , 

-ni  —  1  —  <  1  —  + 


«H2,eff,z 


-1 


^eADLMACHcACH  /  CXP  lvpZ^ff/SADL  j 


«H2,eff,z 


(34) 


On  substitution  for  o  f  ^ach-dl  from  Eq.  (34)  into  Eq.  (32),  the  following  equation  was  obtained  to  predict  o>\  [2  as  a  function  of  the 
distance  z  in  the  anode  reaction  layer: 


(Wh2  =  ot 


cosh 


'FaZ 


t,  sinh 


'FAz 


(35a) 


which  is  valid  for  <5ACH  DL  <  z  <  <5ACH  DL  RL  and  where 
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and 


{1  -  «A™  +  nH2  effz /  kc,H2  exp (m2,eSJpD^/S^) 


„ADL fl^ACH „ACH 


,ADL  /oADLa 


cosh 


-DL 


1  —  tanh 


pA ^ACH— DL 


tanh 


^a^ach-dl-rl 


(35c) 


where  nH2,eff,z  =  igeomMH2/2F. 

The  prediction  of  co\  \2  as  a  function  of  z  at  an  x-plane  requires  the  data  on  «H2,eff,j  corresponding  to  a  geometric  current  density, 
transport  and  electrode  electrochemical  kinetic  parameters,  and  dimensions  of  the  anode  flow  channel  and  diffusion  and  reaction 
layers  parallel  to  the  z  direction.  Also,  the  effective  diffusivity  of  hydrogen  in  the  anode  reaction  layer,  Wfj^eir,  appearing  in  Eq. 
(31b)  is  given  below: 


D\ 


ARL 
H2,eff  : 


t^ARL  c  ARL 

un2  £g 

tarl 
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Eq.  (35a)  was  used  to  obtain  (— d&>R2/dz|.=,5ACH-DL)  and  the  result  substituted  into  Eq.  (27)  to  obtain: 
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(37) 


If  at  an  x-plane,  the  required  geometric  current  density  is  zgeom,  the  hydrogen  mass  flux  that  must  enter  the  anode  reaction  layer  at 
-  _  ^ACH-DL  _  ^ACH  +  ^ADL  js 


A  ARL 
nH2,eff,z 


- 

Jgeom  ^ 
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On  equating  the  right-hand  sides  of  Eqs.  (37)  and  (38),  one  obtains  the  following  equation: 


fi=aJFA 
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/  ^A^ACH-DL 


sinh 
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One  should  solve  Eq.  (39)  for  y  FA  for  a  given  value  of  igeom  using  a  method:  trial  and  error,  Newton’s  or  graphical  method;  or  any 

commercial  computer  software.  One  can  then  determine  the  anode  reaction  overvoltage  \)]A\  from  Eq.  (31b)  corresponding  to  the 
geometric  current  density,  z„eom  ■ 

Chemical  species  mole  balances  in  the  anode  channel  flow  are  now  considered.  Under  the  set  of  operational  conditions  of  a 
reformer,  when  negligible  amounts  of  carbon  monoxide  and  oxygen  are  present  in  the  fuel  gas  supplied  to  a  PEMFC,  the  fuel  feed 
gas  mixture  to  the  cell  anode-side  flow  channel  is  assumed  to  contain  hydrogen,  nitrogen,  carbon  dioxide  and  water  vapor.  It  is 
further  assumed  that  leakage  of  gases  through  the  polymer  electrolyte  separator  from  the  cell  anode  to  cathode  is  negligibly  small. 
The  molar  flow  rates  (mol  s“ 1 )  of  the  chemical  species,  nitrogen,  carbon  dioxide,  hydrogen,  and  water  vapor  at  any  x-location  are 
given  as  follows: 
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(41c) 


for  the  assumption  of 


total  cell  current 


geom 


*geom,ave 


total  geometric  area  of  the  cell  anode 


»,ACI[  _  «,aci I 
iVH20  —  iVH2O,0 


(42a) 


if  the  cell  operational  conditions  are  such  that  the  net  loss  of  water  from  the  gas  mixture  flow  in  the  cell  anode  channel  is  negligible. 
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If  the  net  amount  of  water  transported  from  the  anode  to  cathode  side  of  the  cell  is  £h20>  moles  of  water  per  mole  of  H+  ions 
that  move  from  the  anode  to  cathode  reaction  layer  of  the  cell;  then,  the  molar  flow  rate  of  water  in  the  cell  anode  channel  at  an 
x-location  is  given  as 


lirACH  _  lOACH 
/vh2o  —  7VH2O,0 


?h2oW 


ACH 


^geom' 


dx 


=  Ni 


\rACH 


fH#ACH^. 


h2o,o 


geom.ave 


(42b) 

(42c) 


for  the  assumption  of  igeom  =  'geom.ave- 

The  species  molar  flow  rate  information  provided  in  Eqs.  (40)-(42)  should  be  used  in  Eqs.  (13c)-(13e).  If  the  fuel  gas  mixture 
feed  to  the  cell  anode  channel  contains  small  amounts  of  carbon  monoxide  and  oxygen  in  addition  to  nitrogen,  carbon  dioxide, 
hydrogen  and  water  vapor;  then,  one  may  use  the  following  equations  for  carbon  monoxide  and  oxygen  at  an  x-location  in  the  anode 
flow  channel: 


lirACH  _  lOACH  ,  ArACH 

-,vco2  —  iVco2,o  +  7VCO,0 


vCO2,0  ‘ 


a/ACH  _  »/ACH 
iV02  —  iVO2,0 


i  A/ACH 
2iVCO,0 


(43) 

(44) 


Eqs.  (43)  and  (44)  were  obtained  with  the  assumption  that  the  small  amount  of  carbon  monoxide  present  in  the  feed  would  be 
converted  to  carbon  dioxide  completely  via  the  CO  oxidation  reaction:  CO  +  (1/2)02  CO2,  that  takes  place  preferentially  over  the 
hydrogen  oxidation  reaction  at  the  platinum  surface  for  temperatures  less  than  200  °C  [35].  It  is  here  suggested  that  the  amount  of 
oxygen  present  in  the  fuel  feed  gas  mixture  should  be  such  that  it  is  sufficient  to  completely  convert  CO  at  the  platinum  electrode 
surface  to  CO2,  particularly  for  high  temperature  phosphoric  acid  PEMFCs,  leaving  residual  oxygen  in  the  anode  channel  flow, 
/Vq™,  to  a  negligible  level  to  prevent  possible  oxidation  of  hydrogen. 


2.2.  The  fuel  cell  cathode-side  formulation 


The  flux  of  O2  in  the  ’-direction  through  the  gas  phase  concentration  boundary  layer  prevailing  at  the  porous  surface  of  the  carbon 
cathode  diffusion  layer  (CDL)  facing  the  cathode  channel  flow  is  given  by 


(— ZVoj.eff.z)  =  *c™£gDL(4£H  -  co2,zcdl) 


The  mass  flux  of  oxygen  in  the  ^-direction  through  the  gas  phase  concentration  boundary  layer  is  given  by 
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RT 


(46c) 

(46d) 

(46e) 


As  with  the  anode  channel,  the  cathode  channel  mass  transfer  coefficient  ™  value  can  be  determined  using  semiempirical  equations 
[36,38].  The  gas  mixture  density  in  the  cathode  channel  is  to  be  calculated  from: 


pCCH 

P  =  — (a  =  02,  N2,  H20  mainly) 

a 


(47) 


The  molar  flux  of  H20(V)  in  the  ^-direction  through  the  gas  phase  concentration  boundary  layer  prevailing  at  the  porous  surface  of 
the  CDL  is  given  by 

AT  ,„CCH  CDL/  CCI  I , 

/VH20,eff,’  =  ^C,H2Osg  (cH20,zgDL  _  CH2Q  ) 


(48) 
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The  mass  flux  of  FbCfv)  in  the  z-direction  through  the  gas  phase  concentration  boundary  layer  is  given  by 

•  <7  M  >CCH  CDL  ii^CCH  CCH/  ,  CCH  , 

«H20.eff,z  =  tVH20,eff,zMH20  =  *C,HOeg  M  c  (ftJH20,z™L  “  ^HjO 1 

From  Eqs.  (46a)  and  (49),  one  obtains: 

CCH _ (~ »Q2.eff.z) 

°24DL  02  £CCH£CDLMCCHcCCH 

, ,  _  ,  CCH  , _ «H2Q,eff,z _ 

WH20,z™l  -  H20  +  ytCCHo£CDLMCCHcCCH 

The  nitrogen  molar  flux  along  the  "-direction  is  assumed  zero.  That  is 

-%2,eff.z  =  0-0 

Also 


(49) 


(50) 

(51) 


(52a) 


«N2,eff,z  =  -%2jeff,-MN2  =  0 


(52b) 


The  general  mass  flux  equation,  which  describes  the  mass  transport  of  a  chemical  species  a  in  the  CDL  by  the  mass  diffusion  and 
convection  processes,  is  given  as 


«a,e£f,z  =  ^  a  )  +  COa5^nfre{f'Z  (a,  j8  —  O2,  H20,N2) 


p 


Under  the  steady-state  condition  at  an  x-location: 


(53) 


«02,eff ,z  (at  any  z  within  the  CDL)  =  «o2,eff,zl,=-CRL  =  - 


lgeom 

4F 


Mq2 


The  water  mass  flux  in  the  z-direction  within  the  CDL  at  any  x  is  given  as 

Lcorn  M 1  [2()  (l  [2oigcom  rV/ f  [ 2 ()  .  tgeom  37  J  |20 

«H20,e£f,z  =  - 2F  H - - - — -  =  (!  +  2(h2o) - 2F 

Using  the  expressions  in  Eqs.  (54)  and  (55),  with  «N2,eff,z  =  0: 

y^»j3,eff,z  =  «02,eff,z  +  «H20,eff,z  +  «N2,eff,z 
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Combining  the  information  in  Eqs.  (53)  and  (56c)  to  obtain: 
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where 
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Eq.  (57a)  is  transformed  to 
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(58) 
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The  analytical  solution  to  the  differential  Eq.  (58),  with  the  boundary  condition  at  z  =  ZfsUL  (far-side  of  CDL  with  respect  to  z  =  0), 
coa  =  a>™L,  is  given  below: 


cm 

°>ct  =  (Oaf,  exp 


(1  +  18£H2o)0geom4^H2/2-E)  .  CDL\ 

- TTIcdl - (Z-Zfs  ) 
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\  PUa,eif 


(59) 


valid  for  a  chemical  species  a  in  the  z-range,  z^DL  —  Z  <  z£DL;  a  =  CE,  EhOfv),  Nt. 

Eq.  (59)  is  applied  to  the  chemical  species,  N2,  O2,  and  H20(v),  respectively,  in  conjunction  with  Eqs.  (52b),  (54)  and  (55),  to 
obtain: 


cni 

<wNz  =  exp 


(1  +  18£H2oX!geom4^H2/2C) ,  CDL, 

- TT^CDL - ) 

^N2,eff 


(60) 


rot 

«02  =  «o2,fs  exP 


(1  +  18fH90)(l'geomMH2/2f)  CDL, 

- TTiCDL - ) 

Pu02,ttt 


,  .  (1  +  18()H20)(igeomiWH2/2C  .  CDL,  I 

i  —  exp  ,,rni  IZ  —  Zfs  ) 


pD\ 


,CDL 

02,eff 


8 


1  +  18Ch2o 


(61) 


CDI 

«H20  =  «H20,fs  eXP 


(1  +  18^H2o)(igeom4^H2/2E)  CDL\ 

- TTIcdl - - - (z  “  zfs  ) 

^H20,eff 


1  +  2Ch2Q  \ 

1  +  18^h2o/ 


1  I  (1  +  18^H2o)0geom-^H2/2E)  _CDL, 

1  exP  j  CDL  'Z  4fs  > 

(  ^^H20,eff 


(62) 


The  mathematical  expressions  for  (— ho2,eff,z)  and  «ilo.cit,-  from  Eqs.  (54)  and  (55),  respectively,  are  substituted  into  Eqs.  (50)  and 
(51)  to  obtain  the  expressions  for  <w0i  „cdl  and  &>H?0  _cdl.  The  resulting  expressions  for  &>0 2  _cdl,  <wH90  -cdl,  and  ojn2  -cdl  = 
are  substituted,  respectively,  into  Eqs.  (61),  (62)  and  (60)  to  obtain  the  following  expressions: 


&>02  = 


,CCH 


(igeomMo2/^F) 


02  ^£CDLMCCHcCCH 


exp 


(1  +  18fH2o)(igeomMH2/2E)/  CDL\ 
- (Z  Zfs  ) 


pD\ 


CDL 

02,eff 


1  +  18^h20 


1  J  (1  +  18£H2o)(tgeom^H2/2 F).  CDL,  I 

1  exP  j  ^cdl  '  Zfs 

l  Pu  02,eff  J 


(63) 


(WH2O  = 


CCH  1  (1  4“  2(| [2o)0gcorn /l. F) 

,CCHo£CDLMCCHcCCH 


wH20 


exp 


(1  +  18fH2o)(*geom4/H2/2T')  CDL\ 

- ^5CDL - (Z  —  Zfs  ) 


H20,eff 


J  (1  +  18£H'>o)0geom4^H2/2-E)  CDL, 

exP  \ - ;,ncm. - (z  -  zfs  ) 

(  PuU20,efi 


1  +  2Ch2Q  \ 

1  +  18^h2o/ 


(64) 


CDI 

"Nz  =  «N2,fS  eXP 


(1  +  18(lH2o)(*geomiWH2/2^)  CDL, 

(z.  Zfs  ) 


PD 


,CDL 

N2,eff 


where 


(65a) 


,  ,CDL  _  ,  CCH 
<wN2,fs  —  ^2  _ 


1  gcom 


4  ^7eCDL4^CCHcCCH 


(2  +  4fH2o)^§? 


Mo2 

lCCH 


VC,H20  C,02 


(65b) 


The  Eqs.  (63)-(65a)  can  be  used  to  predict  the  mass  fraction  profiles  of  oxygen,  water  vapor  and  nitrogen  in  the  gas-filled  pores  of 
the  cathode  diffusion  layer  (CDL)  for  z„sDL  <  z  <  z^DL-  The  mass  fractions  of  the  chemical  species  oxygen,  water,  and  nitrogen  at 
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the  interface  between  the  CRL  and  CDL,  but  on  the  CDL-side,  are  obtained  by  setting  z  =  ZnSDL,  with  (T™1*  —  z£DL)  =  —  <5CDL. 


wo2.z£PL 


CCH 

^o2 


( igeoin  M  O2  /  A  F) 

kCWECDLMCCHcCCH 


exp 


(1  +  18?H2o)0geom^H2/2.F)  „CDL 

oOCDL  d 

PL>02,eff 


8 


1  +  18£h2o 


J  “(I  +  18^H2o)(igeom^H2/2f)  CDL 

”eXPl  P°  &  j 


(66) 


h2o^ 


CCH  (1  "3"  2^| [2())figcoin 41} \2q/2F) 


uH20 


^c™oegDLMCCHcCCH 


exp 


"(1  ~t~  18£H2o)0geom^H2/2F)  ^DL 

/,/)CDI  6 

Pu  H20,eff 


1  +  2^h2Q 
1  +  18^h2o 


1  J  (1  +  18fH->o)((geom4^H2/2F)  SCDL 

XP  I 

L  Pu  H20,eff 


(67) 


<iZcDL=4°fsexp 


-(1  +  18CH2o)(JgeomMH2/2F)5CDL' 


n  nCDL 
P^N2,eff 


(68) 


In  the  cathode  reaction  layer,  oxygen  is  transported  as  well  as  reduced  via  the  electrochemical  reaction  given  by  Eq.  (2)  at  the  catalyst 
platinum-electrolyte  active  sites-(s)  to  produce  water.  If  the  first  electron  transfer  step,  CF  —  s  +  H+  +  e_  HO2  —  s,  is  assumed  to 
be  the  rate-limiting  step  in  the  kinetic  mechanism  of  the  overall  oxygen  reduction  process,  the  oxygen  reduction  rate  per  unit  active 
surface  area  of  the  electrocatalyst  is  given  by 


(— ro2)  = 


—  J  exp 
4  F  V  RT 


(69) 


where  the  cathodic  exchange  current  density,  /q  s,  is  given  by 


*0.s  =  «0,s,7b,po2,0  eXP 


R  V  Tq 


po2 

Po2,o 


(70) 


For  the  first  electron  transfer  step  as  the  kinetic  rate-limiting  elementary  step,  the  reaction  order,  n  =  I  [41]: 


\rf\  =  I EC 


Ec  I  — 

-^revl 


(A </>)C  -  (A </>rev)C 


—  I  (0Pt  (/’electrolyte)  (<APt  •/’electrolyte  )ri 


(71) 


where  hS  is  the  actual  voltage  at  the  cathode  reaction  layer  (V). 

For  the  assumption  of  spherical  crystallites  of  the  platinum  catalyst,  with  average  size,  /  pt  (radius);  the  platinum  surface  area  per 
unit  mass  of  the  catalyst  is  given  by 


sc  = 


'PtPPt 


(72) 


If  some  fraction  of  the  platinum  catalyst  surface  is  not  in  contact  with  the  electroactive  regions  of  the  polymer  electrolyte;  then,  the 
reaction  effective  electrocatalyst  interfacial  area  is  given  by 


X  _  c  c 

5eff  —  aetfs 


(73) 


If  the  catalyst  Pt  loading  in  the  cathode  reaction  layer  is  nipt,  the  electrochemical  reaction  active  electrocatalyst  surface  area  per  unit 
cathode  reaction  layer  volume  is 


c  c  c 
a  =  ,seffWpt 


(74) 


The  relation  of  /«pt  to  the  catalyst  loading  per  unit  geometric  area,  Wpt  geom,  of  the  cathode  reaction  layer  is  given  as 


tPt,geom 


=  sCmFv  [(g  Pt)cmge2om] 


(75) 
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The  molar  oxygen  consumption  rate,  via  the  cathodic  electrochemical  reaction,  per  unit  cathode  reaction  layer  (CRL)  volume  is 
given  by 


c 

a  exp 


(76) 


The  oxygen  consumption  rate,  in  mass  units,  in  the  CRL  is  given  by 


dC 

^02, mass 


ac  Mq2  exp 


(77) 


Oxygen  diffuses  toward  the  solid  polymer  electrolyte  separator  layer  along  z-coordinate  in  the  direction  of  decreasing  z  while  being 
reduced  by  the  electrochemical  reaction  at  the  electrocatalyst  interface  effective  active  sites.  Based  on  the  coupled  oxygen  transport 
and  electrochemical  reduction  processes,  a  differential  model  equation  describing  the  oxygen  mass  fraction  variation  with  distance 
Z,  obtained  by  the  application  of  steady-state  oxygen  mass  balance  over  a  spatial  element,  is  given  by 


(-R 


c 

(>2,mass 


) 


(78) 


In  the  derivation  of  Eq.  (78),  the  product  of  gas  mixture  density  and  effective  oxygen  mass  diffusivity,  i.e.  pOo^eff’  was  assumed 
invariant  with  respect  to  z  in  the  thin  cathode  reaction  layer  (CRL).  The  gas  mixture  density,  p,  may  be  taken  equal  to  that  in  the 
cathode  channel  gas  mixture  at  an  x-location  as  a  first  approximation.  The  gas  mixture  density  at  an  .r-location  can  be  computed 
from  Eq.  (47).  The  effective  oxygen  mass  diffusivity,  in  the  cathode  reaction  layer  is  given  by 


D 


CRL  _ 
02,eff  — 


D 


CRL  CRL 
02(g)£g 
_CRL 
rp(g) 


(79) 


and 


CRL  _  ,  /  CRL  ,  CRL  ,  CRL  ,  CRLx 

fcg  —  1  ffcpt  T  c-q  T  fcSPE  T  other! 


(80) 


If  the  gas-filled  average  pore  radius  in  the  CRL  is  less  than  100  A,  Dq^  should  be  corrected  for  the  Knudsen  effect. 

Assuming  a  first-order  reaction  with  respect  to  oxygen  (n  =  1)  and  using  the  expression  for  io,s  in  Eq.  (70);  the  oxygen  consumption 
rate  per  unit  CRL  volume,  Eq.  (77)  leads  to 


(-*! 


02,mass 


)  = 


'  io,s,TQ,po2,oMCCHac 

4 F 


exp 


zT_ 

R 


1 

% 


■C  i  „C  i  p 


RT 


a>o2 


Substitution  for  (—  Rq  mass)  from  (81)  into  (78)  leads  to 


d2a>p, 

dz2 


Ffa>o2 


where 


CCH 


n  = 


(%s,T0,PO2aCM' 

{  4F(pCCH^CRLff) 


exp 


ZSL 

R 


1 

% 


|  rf  |  F 
RT 


(81) 


(82a) 


(82b) 


For  the  isothermal  condition,  with  the  assumption  of  constant  |);c|  in  the  relatively  thin  cathode  reaction  layer,  Eq.  (82a)  was  solved 
using  the  following  boundary  conditions: 


CRL 


and 


z  =  z 


,  _  „CRL 
Z  Zfs  , 


(  d&>o2 


V  dz 


=  0 


, ,  _  ,  CRL— 

a>o2  —  c°o2,fs 


(83a) 


(83b) 


cnt  _ 

where  a>(h  fs  =  mass  fraction  of  oxygen  on  the  negative  side  of  the  interface  between  the  CRL  and  CDL.  The  solution  to  the 
differential  Eq.  (82a)  is  given  as 


&>02  =  oq 


cosh 


Fj  z  )  -  Ci  sinh  W  F,  z 


(84a) 
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where 


co, 


ai  = 


CRL— 

02.fs 


cosh 


PC,CRL 
r  1  <-fs 


1  —  tanh  I  J  FpzpRL  )  tanh 


pC_CRL 
M  '■fs 


(84b) 


valid  for  zpsRL  <  z  <  ZpRL. 

In  general  (activity  of  oxygen  on  the  negative-side  of  the  interface  between  the  CRL  and  CDL)  =  (activity  of  oxygen  on  the 
positive-side  of  the  interface  between  the  CRL  and  CDL),  with  the  assumption  of  equilibrium  with  respect  to  oxygen  at  the  interface. 
However,  for  the  presence  of  gas  phase  in  the  pores  of  CRL  and  CDL,  it  is  quite  reasonable  to  assume: 


CRL— 

0O2,fs 


(85) 


Using  the  information  provided  for  cot  CDL  in  Eq.  (66)  in  conjunction  with  Eq.  (85),  Eq.  (84b)  becomes: 

02 » 2ns 


a\  = 


-  (igeomMo2/4F)/^CH£CDLMCCHcCCH]exp[(1  +  18((H2o)0geomMH2/2F),5CDL / ptig£B\ 


cosh 


rC.CRL 
M  zfs 


1  -  tanh  (  J  FpzpsRL  )  tanh  (  J  FpzpRL 


(8/(1  +  18?h2o))[1  -  exp{-(l  +  18CH2o)0geomMH2/2F)3CDL/(pD^ff)}] 


->CDL 


cosh 


pC  CRL 
M  '■fs 


1  —  tanh 


*fz2*L)  tanh  (^f4RL) 


(86) 


=  tanh 


rC  CRL 
rl  '■ns 


(87) 


Eq.  (84a)  can  be  employed  to  predict  the  oxygen  mass  fraction  profile  in  the  cathode  reaction  layer  as  a  function  of  distance  z.  At 
an  x-plane,  the  oxygen  mass  flux  along  the  z-direction  is  given  by 


(  ■  ,  _  CCH  nCRL 

(  nC>2,eff,z)  —  P  D, 


dtt>o2 


02,eff 


dz 


(88a) 


tgeom 


Mq2 


4  F 


(88b) 


Using  the  information  provided  for  ao2  in  Eqs.  (84a),  in  conjunction  with  Eqs.  (88a)  and  (88b),  the  following  result  is  obtained: 


^geom 


Mq2 


=  p 


CCH  nCRL 

Ur 


4/Z  02,eff 

Eq.  (89)  can  be  re-arranged  as 


a  l  y  7*1 


sinh 


if  Z -  Ci  cosh 


(89) 


tgeom 


Mq2 


4FpCCHDCRL^ 


O'  1 


sinh 


^1C4RL  -Cl  cosh  K/ifzgRL 


(90) 


For  the  steady  operation  of  a  fuel  cell,  Eq.  (90)  can  be  solved  for  by  a  trial-error  or  graphical  procedure  for  a  fixed  value  of  LCOm 
at  an  x-plane.  Then,  using  the  defining  Eq.  (82b)  for  Fp,  one  can  obtain  the  cathodic  electrochemical  reaction  activation  voltage 
loss,  |i7C|,  for  a  fixed  igeom  value  at  an  x-plane.  Also,  one  can  find  the  average  cathodic  activation  voltage  loss,  |i;c|ave  using  the  cell 
average  geometric  current  density.  This  rigorous  mathematical  formulation  shows  how  one  can  obtain  the  activation  voltage  loss  at 
the  cathode  reaction  layer  by  taking  into  account  the  effect  of  mass  transport  of  oxygen  through  the  various  layers  of  the  cell  on  the 
CRL  side. 

The  chemical  species  mole  flow  rates  at  the  cathode  channel  inlet,  x  =  0  plane,  are  given  by  /Vp(pH,  (i  =  CL,  H2O,  N2).  Under  the 
cell  steady-state,  isothermal  operational  conditions,  the  chemical  species  molar  flow  rates  in  the  cathode  channel  at  an  x-plane  are 
given  as  follows: 


yCCH  _  iha.ii 
v  —  %n 


v02,x 


<rCCH 
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CCH 


o2,o 


4  F 


fcgeom 


dx 


(91a) 
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with  igeom  —  igeom, ave- 


X/CCH  _  a  rCCH 

iVo2,jf  —  yvo2,o 


X/CCII  _  a rCCH 
/VH20,x  —  /VH2O,0 


W 


CCH , 


cgeom,ave 


4  F 


-  Eh2o  \ 


2  F  J 


W 


CCH 


tgeom 


dx 


—  +  (1  +  2^H2o) 


X/CCH  _  a2CCH 
7VN2  —  a'n2,o 


lyCCI  ^  x/oaongaw 

2 F 


The  total  gas  mixture  molar  flow  rate  in  the  cathode  flow  channel  at  an  x-plane  is  given  by 


OCCH 


n: 


=  E^CH  =  E^CH  +  (1  +  4^h2o)- 


/CCH  Jo  ;geom  dx 


4F 


and,  with  Igeom  —  *geom,ave  • 

=  EE?H  +  0  +  4Ch2o) 


W 


CCH. 


‘geom.ave 


4F 


Mole  and  weight  (or  mass)  fractions  of  chemical  species  at  an  x-plane  in  the  cathode  flow  channel  are,  respectively,  given  by 
aT'CII 

CCH  _  lyi.x 
yi  /V;CC" 

OCCH 


N 


l,X 


E/Eo  +  (VTCCH/4F)(1  +  4Ch2o)  fj  igeom  dx 

with  igeom  =  igeom, ave  3nd  i  =  O2,  HgOfy),  N2 : 


yi 


CCH 


frCCH 

i,x 


E/Eo  +  (1  +  4CH2o)(WCCH*igeom,ave/4F) 
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EiN^Mi 
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l,X  1 


Ei*TMi  +  WCCH(1  +  ?H2oMh2o)  /o  igeom  dx/F 

with  igeom  =  igeom.ave  Und  i  =  O2,  hEOfy),  N2 ! 


,CCH 


^ccuM 

T  1.X  l 


E  Mr  Mi  +  (1  +  ?H2oMH2o)(WCCH^igeom,ave/i7) 


(91b) 

(92a) 

(92b) 

(93) 

(94a) 

(94b) 

(95a) 

(95b) 

(95c) 

(96a) 

(96b) 

(96c) 


The  oxygen  weight  fraction  at  an  x-plane  in  the  cathode  flow  channel,  as  provided  by  Eq.  (96b)  or  (96c),  is  to  be  used  in  Eq.  (86)  to 
calculate  the  parameter  a  \  for  its  subsequent  use  in  Eq.  (90)  to  calculate  the  parameter  I:\  ;  hence,  the  activation  voltage  loss,  if  at 


an  x-plane.  Eq.  (96c)  can  be  used  to  compute  the  oxygen  mass  (or  weight)  fraction,  &>qCiJ,  as  a  function  of  distance  x  in  the  cathode 


it’ 


o2,*> 


channel.  Then,  the  average  weight  fraction  of  oxygen,  of\  ave,  can  be  calculated  from  the  following  equation: 

r  rrw 
rx=L 

Jx= 0 


rx=Lccil  CCH 


CCH 
^02, ave 


«5r dA 


lcch 


(97) 


The  average  oxygen  weight  fraction,  ®o™ve,  should  then  be  used  to  calculate  the  average  activation  voltage  loss  at  the  cathode 
reaction  layer,  |  ?;c  | ave ,  for  an  average  geometric  current  density  of  igeom,ave  using  Eqs.  ((90),  with  igeom  =  igeom.ave),  ((86),  with 


off  —  ®0?,ave’  igeom  —  igeom.ave);  Und  ((82b),  with  | T)  \  —  1)  lave)- 
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3.  Application  of  the  developed  formulation 

The  formulation  provided  in  Section  2  may  be  employed  to  calculate  an  average  cell  voltage  at  a  given  or  desired  average  current 
density  for  a  hydrogen/air  PEMFC,  using  an  acid  electrolyte  medium,  for  example,  Nafion-Teflon-Zr(HP04)2  composite [41], 
PBO/PBI  doped  with  phosphoric  acid,  or  sulfonated  poly(arylene  ether  sulfone)  [42].  Example  cell  operational  conditions  can  be  an 
elevated  temperature  of  1 10  °C,  1  atm  total  pressure,  with  the  water  vapor  partial  pressure  less  than  saturated  vapor  pressure  in  the 
anode-side  hydrogen  fuel  stream  and  50%  or  even  lower  relative  humidity  in  the  cathode-side  oxidant  air  stream  to  maintain  oxygen 
partial  pressure  at  a  sufficiently  high  level  to  prevent  excessive  activation  voltage  loss  at  the  cathode  reaction  layer.  One  should 
determine  the  total  ohmic  voltage  loss,  ifn  ,  for  a  desired  average  current  density,  igeom.ave,  using  Eqs.  (7)— (12a).  The  average 
activation  voltage  losses,  |  lave  and  b)C|ave  can  be  calculated  using  Eqs.  (39)  and  (90).  The  average  reversible  cell  voltage,  E*™?ave, 
should  be  calculated  using  Eq.  (4)  at  the  arithmetic  average  of  the  partial  pressure  values  of  hydrogen,  oxygen,  and  water  vapor  at 
the  anode  and  cathode  channel  inlets  and  exits.  Then,  the  cell  actual  average  voltage  can  be  computed  from  the  following  equation: 

C“  =  £revUave  “  4,ave  ~  (l/lave  +  lbC|ave)  -  if*  (98) 

It  is  here  noted  that  the  mass  transport  related  voltage  loss,  rfn,  should  be  set  equal  to  zero  if  the  E“".lve  is  computed  at  the  arithmetic 
average  of  the  partial  pressure  values  of  hydrogen,  oxygen,  and  water  vapor  at  the  anode  and  cathode  channel  inlets  and  exits  as 
they  prevail  at  the  anode  and  cathode  reaction  layers.  If  £“1* ave 's  computed  at  the  average  of  these  species  partial  pressures  as  they 
prevail  in  the  channel  flow  streams;  then  there  may  still  be  some  need  to  apply  a  mass  transport  correction  through  the  use  of  an 
empirical  equation  such  as  [39]: 

r)m  =  m  exp(m'geom)  (99) 

where  m  =  2. 1 1  x  10-5  V,  n  =  8  cmj;eom  A- 1 ,  ige om  in  (A  cm^eom). 

The  effect  of  hydrogen  fuel  crossover  from  the  cell  anode  to  cathode  on  the  cell  open  circuit  voltage  is  significant  [39].  For 
example,  if  the  hydrogen  crossover  flux  corresponds  to  the  internal  current  density  (=(no.  of  g  mol  of  hydrogen  diffusing  from  the 
anode  to  cathode  reaction  layer,  through  the  electrolyte  separator  layer,  per  unit  geometric  area  per  s)  x  (2 F  coulomb  per  mole  of 
hydrogen  crossover))  =  1  (mA  cm;eom),  the  open  circuit  cell  voltage  is  about  0.97  V  instead  of  the  thermodynamic  reversible  voltage 
of  1.2  V  at  25  °C.  In  the  practical  current  density  range  of  igeom  >  100mAcm“eonl,  the  effect  of  hydrogen  crossover  on  the  cell 
operating  efficiency  is  not  of  great  significance  depending  upon  the  diffusion  characteristics  of  the  membrane  separator  layer.  One 
can  compute  the  voltage  loss  due  to  the  hydrogen  crossover  flux  using  the  simple  approach  given  in  Ref.  [39]. 


4.  Single  cell  electric  power  and  thermal  efficiency 

Corresponding  to  a  desired  average  geometric  current  density,  Zgeom,ave>  the  total  cell  current  (A),  is  given  by 

rcell  _  •  /  j  ACHti/ACH\ 

—  *geom,ave  vv  ) 


(100) 


where  LACH  and  WACH  are  the  dimensions  of  the  anode  electrode  parallel  to  the  x-  and  y-coordinates  as  shown  in  Fig.  1.  The 
magnitude  of  the  cell  electric  power  (W),  is  given  by 


I  Wc,e11  .  I  =  ycenrc 
•  electric  l  ave  t 


cell  rcell 


Hydrogen  consumption  rate  in  the  cell  (mol  s  *): 


rcell 


(—A  AT,6,11)  = 

v  h2  >  2F 

Rate  of  thermal  energy  production  if  the  entire  amount  of  hydrogen  consumed  in  the  cell  were  combusted: 


ARtceU  =  (-AN^l)AHT  = 


rcell 

b 


A  Ht 


(101) 


(102) 


(103) 


where  A Hj  is  the  enthalpy  of  the  hydrogen  oxidation  reaction,  Eq.  (3),  at  the  cell  temperature,  T  (K).  For  the  gas  mixture  ideal 
gas  behavior,  A  Ht  —  A  HT  =  Yi  v,-  / (  V  —  standard-state  enthalpy  of  the  reaction,  Eq.  (3),  at  the  cell  temperature,  T,  per  gmol  of 
hydrogen  oxidized  to  water  vapor  (J  mol-1).  Eq.  (103)  now  becomes: 


S.S.  Sandhu,  J.P  Fe line r  /  Journal  of  Power  Sources  161  (2006)  1133-1153 


1151 


AH}(]mo\  *)  is  given  by  [36]: 


AH°t  =  5>A  HliJo  +  R 


A  B  7 

AMT-  7o)  +  —  (T2 


'f  A  C  O  Q 
A?)  +  -^(r3  -  r3) 


,  1  l 

AD - 

T  T0 


(105) 


where  R  =  8.314  (Jmol_1  K-1);  7o  =  298.15 K;  vH2(g)  =  -1,  vo2(g)  =  -1/2,  yH20(g)  =  1;  &hIh2,T{)  =  A hio2.t0  =  0  (Jm°l~  )> 
A77f°H2o,Tb  =  —  24 1818  (Jmor1);  AA  =  X>,-A/  =  -1-5985,  A B  =  =  0.775  x  lO^KT1,  AC  =  X>;C;  =  0-0  K“2, 

AD  =  vi  Oi  =  0.130421  x  105  K2,  where  A/,  fi,,  C,,  D,  are  the  heat  capacity  coefficients  of  an  ideal  gas  i. 

The  single  cell  thermal  efficiency  is  defined  as 


cell  _ 

'/thermal 


- 1  wc,dl  - 

1  electric  I 
z>cell 


2  FV! 


cell 


A/7tce"  {-AH}) 

Thermal  energy  that  must  be  removed  from  the  cell: 

I  2 remove  I  =  d  ~  'SmalX- A77tcdl) 

to  maintain  the  cell  at  a  desired  temperature,  T. 


(106) 


(107) 


5.  Extension  of  the  single  cell  formulation  for  the  performance  prediction  of  a  stack  of  identical  cells  connected  in  series 
(steady-state,  isothermal  operation) 

It  is  here  assumed  that  the  anode-side  flow  channels  of  the  individual  PEMFCs  connected  in  series  are  fed  with  the  identical  rates 
of  fuel  (hydrogen  or  a  gas  mixture,  containing  hydrogen,  from  a  fuel  reformer).  The  cathode-side  flow  channels  of  the  individual 
PEMFCs  are  also  fed  with  the  identical  rates  of  oxidant,  air.  The  fuel  and  oxidant  (air)  flows  are  in  the  same  (i.e.  parallel)  flow 
direction.  For  the  cells  connected  in  series,  the  same  amount  of  total  current  passes  through  each  cell.  That  is  to  say  that 


rcell  _  rStack 

n  —  yt 


(108) 


where  /(stack  —  total  current  from  a  stack  of  PEMFCs  connected  in  series.  For  n  PEMFCs  connected  in  series  in  a  stack,  the  total 
stack  voltage  is  given  by 


T/Stack  cell 

vt  n  *ave 


(109) 


where  a  single  cell  average  voltage  is  obtained  from  Eq.  (98).  The  cell  stack  electric  power  output  is  given  by 


I  Electric 


T /Stack  /Stack  T/cell  /Stack 
=  yt  A  =wVaveA 


(110) 


The  hydrogen  fuel  consumption  rate  per  cell  which  is  equal  to  the  rate  of  decrease  in  hydrogen  molar  flow  rate  in  the  cell  anode-side 
flow  channel  is  given  by 


-AiVH2  = 


rcell 


2  F 


/stack 
ft _ 

2  F 


(111) 


The  total  hydrogen  fuel  consumption  rate  in  the  cell  stack  of  n  cells  connected  in  series  is  given  by 


-AAA2  =  «(-A!Vh2)  = 


nl\ 


stack 


2  F 


(112) 


The  thermal  energy  production  rate  is  given  by  the  following  equation  if  the  entire  amount  of  hydrogen  fuel  consumed  in  the  cell 
stack  were  combusted  (i.e.  oxidized  to  water  vapor): 


A77tstack  =  (— AIVjj2)A/7y 


nl. 


stack 


2  F 


AH°T 


(113a) 

(113b) 


The  cell  stack  thermal  efficiency  is  given  by 


stack 

''thermal 


_  I  ti/stack  I 
I  electric  I 
Attack 


2  FVl 


cell 


{-AH}) 


^thermal 


(114) 


1152 


S.S.  Sandhu,  J.P.  Fellner  /  Journal  of  Power  Sources  161  (2006)  1133-1153 


Notice  that  the  cell  stack  thermal  efficiency  is  equal  to  a  single  cell  thermal  efficiency  for  a  set  of  identical  single  cells  connected  in 
series  operating  under  the  isothermal,  steady-state  conditions.  Total  heat  removal  rate  from  the  cell  stack  to  maintain  its  temperature 
at  a  desired  level  of  T  is  given  by 


=  (1  -  '£mal)(-A//tStaCk) 


(115a) 


=  (1  — CCrLt) 


n/tstack 

2  F 


(-A  H°) 


(115b) 


where  n  is  the  number  of  identical  PEMFCs  connected  in  series. 

The  total  molar  flow  rate  of  a  chemical  species  i  at  the  cell  stack  inlet  to  provide  fuel  to  the  anode-side  channel  of  each  cell  in 
the  stack  is  represented  by  A^H,stack  (i  =  H2,  PEO^)  for  humid  hydrogen  feed  or  i  =  H2,  PEO^),  CO2,  N2,  etc.  for  fuel  feed  to  the 
cell  stack  from  a  fuel  reformer).  Total  molar  flow  rate  of  a  chemical  species  i  at  the  cell  stack  inlet  to  provide  oxidant  (air)  to  the 
cathode-side  channel  of  each  cell  in  the  stack  is  represented  by  N^n  ,  ,stac  (i  =  O2,  PEO(V),  N2  for  humid  air  feed).  Total  chemical 
species  molar  flow  rates  at  the  cell  stack  fuel-  and  oxidant-side  exits  are  given  as  follows: 


ACH.stack  n  A 


stack 


irACH^tsck  -jCtj 

^Vexit,H2  ^'in,H2  2^F~ 


»rACH,stack  VrACH, stack  _ 

^exit,H20(v)  ^'in,H20  n 


rstack 


Ch2o 


^ACH,stack  =  ^ACH,stack  (-  =  ^  etc  ) 

stack 


VtCCH. stack  _  VrCCH.stack  _  nh 
^exit,02  ^in,02  4/7 


»rCCH, stack  _  VrCCH, stack  , 

■'vexit.H20(v)  —  JVin,H20(v)  '  n 


VrCCH, stack  VrCCH, stack 

iVexit,N2  ~  7Vin,N2 


(1  +  2feo) 


rstack 

M 

2  F 


The  total  gas  mixture  molar  flow  rates  at  the  stack  channel  exits  are  given  as  follows: 


N. 


ACH,  stack 


exit.t 


EVrACH, stack  _  Vr AC 

^  exit, i  /  J'in.i 


ACH,  stack 


(1  +  2Ch,o)-; 


rstack 


2  F 


N. 


CCH,  stack 


exit,t 


EVtCCH, stack  _  v, CC 
^exit ,i  /  j'in.i 


CCH,stack 


(1  +  4^h2o)A 


rstack 


AF 


( i  —  H2,  HaOiv),  CO2,  N2,  etc.) 

(/'  =  O2,  H20(V),  N2) 


(116a) 

(116b) 

(116c) 

(116d) 

(1 16e) 
(1 16f) 

(117a) 

(117b) 


6.  Concluding  remarks 


Development  of  a  mathematical  model  to  simulate  the  performance  of  a  PEMFC  at  an  elevated  temperature  has  been  presented. 
The  presented  model  is  deemed  sufficiently  rigorous,  yet  not  unwieldy  for  its  application.  Effect  of  the  reactant  species  mass 
transport  on  the  electrochemical  kinetic  reaction  rates  in  the  cell  anode  and  cathode  reaction  layers  have  been  properly  accounted  for 
to  correctly  predict  the  electrochemical  reaction  overvoltages  or  activation  voltage  losses  in  the  cell  thin  reaction  layers.  The  single 
cell  developed  model  formulation  has  been  linked  to  the  stack  model  equations  to  predict  the  performance  of  a  stack  of  a  number 
of  identical  PEMFCs  connected  in  series.  Provided  the  sufficiently  accurate  values  of  the  mass  transport,  charge  transport,  and 
electrochemical  kinetic  parameters  are  made  available,  the  presented  model  formulation  can  be  employed  to  numerically  simulate 
the  performance  of  a  single  PEMFC  as  well  as  of  a  stack  of  a  number  of  identical  PEMFCs  connected  in  series,  utilizing  an  electrolyte 
such  as  sulfonated  polyarylene  ether,  PBO/PBI  doped  with  phosphoric  acid,  orNafion-Teflon-Zr(HP04)2,  at  an  elevated  temperature 
(>100°C)  for  the  isothermal  and  steady-state  operating  conditions. 
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